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Solutions obtained by the quasilinearization method (QLM) are compared with the WKB solu- 
tions. Expansion of the p-th QLM iterate in powers of h reproduces the structure of the WKB 
series generating an infinite number of the WKB terms with the first 2'' terms reproduced exactly. 
The QLM quantization condition leads to exact energies for the Poschl- Teller, Hulthen, Hylleraas, 
Morse, Eckart potentials etc. For other, more complicated potentials the first QLM iterate, given 
by the closed analytic expression, is extremely accurate. The iterates converge very fast. The sixth 
iterate of the energy for the anharmonic oscillator and for the two-body Coulomb Dirac equation 
has an accuracy of 20 significant figures. 



PACS numbers: 03.65.Ca, 03.65.Ge, 03.65.Sq 



I. INTRODUCTION 



The quasilinearization method (QLM) and its iter- 
ations were constructed nl as a generalization of the 
Newton-Raphson method [J for the nonlinear differen- 
tial equations to yield rapid quadratic and often mono- 
tonic convergence to the exact solution. It does not 
rely on the existence of any kind of smallness parame- 
ter. The derivation of the WKB solution starts by cast- 
ing the Schrodinger equation into nonlinear Riccati form 
and solving that equation by expansion in powers of h. 
It is interesting instead to solve this nonlinear equation 
with the help of the quasilinearization method (QLM) 
whose application to physical problems are discussed in 
works HtHilMQ and compare with the WKB results. 

The goal of this work is to point out that QLM iter- 
ates which are expressible in a closed integral form pro- 
vide better approximation than the usual WKB. We show 
that the p-th QLM iterate when expanded in powers of 
H reproduces the structure of the WKB series generating 
an infinite number of the WKB terms with the first 2^ 
terms reproduced exactly. A similar number of the next 
terms are reproduced with approximately correct coeffi- 
cients. We prove also that the exact quantization con- 
dition in any QLM approximation, including the first, 
leads to exact energies for many known physical poten- 
tials used in molecular and nuclear physics such as the 
Poschl- Teller, Hulthen, Hylleraas, Morse, Eckart etc. In 
the general case of arbitrary potentials that do not have 
a simple analytic structure, we illustrate that both the 
wave functions and energies are very well reproduced by 
the first QLM iterate and show significant improvement 
over those obtained by the usual WKB approximation. If 
the initial QLM guess is properly chosen, the wave func- 
tion in the first QLM iteration, unlike the WKB wave 
function, is free of unphysical turning point singularities. 
Since the first QLM iteration is given by an analytic ex- 
pression [3j S S 13 ) it allows one to analytically estimate 
the role of different parameters and the influence of their 
variation on different characteristics of a quantum sys- 
tem. The next iterates display very fast quadratic con- 



vergence so that accuracy of energies obtained after a 
few iterations is extremely high, reaching up to 20 sig- 
nificant figures for the sixth iterate as we show on the 
examples of the anharmonic oscillator and the two-body 
Dirac equation with the Coulomb potential. 

The usual WKB substitution x(r) — 

Cexp(A/ y{r')dr'^ converts the Schrodinger equa- 
tion to nonlinear Riccati form 



dy{z) 



+ {k'{z)+y\z))^Q, 



(1) 



dz 

where k'^{z) = E-V - l{l + l)/z^, A^ = 2™/^^ z = \r. 

The proper bound state boundary condition for poten- 
tials falling off at z ~ zo ~ oo is 2/(2:) = const at z > zo. 
This means that y'{zo) = 0, so that Eq. f^l at z ~ zp 
reduces to fc(zo)^ + j/^(zo)) = or y{zo)) = ±jfc(zo). We 
choose here to define the boundary condition with the 
plus sign, so that y{zo) = ikizn). 

The quasilinearization 0) S Q ^^ this equation gives a 
set of recurrence differential equations 



dypjz) 
dz 



2/p-iW-2yp(z)j/p-i(z)-/c2(z). (2) 



with the boundary condition yp{zo) = ifc(zo). 

The analytic solution .6] of these equations expresses 
the p-th iterate j/p(z) in terms of the previous iterate: 



yp{z) = /p_i(z)- 



, dfp-i{s) 
ds — ^ exp 



fp-i{z) 



ds 

yl_,{z)~k\z) 



Vp~i{t)dt 



(3) 



2yp-i{z) 

Indeed, differentiation of both parts of Eq. © leads im- 
mediately to Eq. ^ which proves that yp{z) is a solution 
of this equation. The boundary condition is obviously 
satisfied automatically. 

The successive integrations by parts of Eq. © lead [g 
to the series 



y,(z) = £/:(f)(z) 



(4) 



C<;r\z) 



with £n'{z) given by recursion relation 
l/(2yp_i(z))d(-4f2,(z))/dz, Ci?{z) = Vi(z). Since 
d/dz — gd/dr, g = A^^ — h/\/2m, Eq. Q represents 
the expansion of the p-th QLM iterate in powers of g, 
that is in powers of h, which one can compare with the 
WKB series. For the zerotli iterate yo{z) it seems natu- 
ral to choose the zeroth WKB approximation, that is to 
set yo{z) = ik{z), which in addition automatically sat- 
isfies the boundary condition. However, one has to be 
aware that this choice has unphysical turning point sin- 
gularities. According to the existence theorem for linear 
differential equations 0, if yp-i{z) in Eq. (0) is a discon- 
tinuous function of z in a certain interval, then yp{z) or 
its derivatives may also be discontinuous functions in this 
interval, so consequently the turning point singularities 
of yo{z) may propagate to the next iterates. 

Eq. © gives an especially simple expression 
[g| for the first QLM iterate yi{z) — ik{z) — 
i J^ dsfc'(s) exp [— 2i/^ k{t)dt\ which thus is expressible 
in a closed integral form. This expression takes into ac- 
count, though approximately, an infinite number of the 
WKB terms corresponding to higher powers of h, as it 
will be shown in the next section. In view of this it is a 
better approximation than the usual WKB. 

To obtain the WKB series one has to expand the so- 
lution y of the Riccati equation QJ in powers of h. This 
is easy to do by looking for y in the form of a series and 
equating terms of identical powers of g, 



y — /_^ y ^m: 



Y„. 






(5) 



The derivatives in this and subsequent expressions are 
in the variable r. The zeroth WKB approximation Yq 
is given by Y^ = ik. The comparison of expansion of 
the first QLM iterate in powers of h and the WKB series 
was originally performed in |y| . There it was shown that 
the expansion reproduces exactly the first two terms and 
also gives correctly the structure of the WKB series up to 
the power g^ considered in these works, generating series 
with proper WKB terms, but with different coefficients. 
Comparison of Eqs. 10} and iQ in the present work shows 
that this conclusion is true also if higher powers of g are 
taken into account. 

The computation of the expansion of the second QLM 
iterate y2 in powers of h is done by reexpanding the term 
1/(2 j/i) in Cn (z) in the series in powers of g and keep- 
ing the powers up to g'^ inclusively in this expression 
as well as in the sum in Eq. 10}. This procedure per- 
formed with the help of Mathcmatica 8] shows that ex- 
pansion of j/2 reproduces exactly already the first four 
terms of the WKB series. It also gives the true structure 
of the next terms of the WKB series, generating series 
with proper WKB terms which have approximately cor- 
rect coefficients. The expansion of j/a is obtained in the 



similar fashion. It reproduces exactly the first eight terms 
of the WKB series, that is all the terms up to the power 
5^" inclusively. 

Summing up, we have proved that the expansion of the 
first, second and third QLM iterates reproduces exactly 
two, four and eight WKB terms respectively. Since the 
zeroth QLM iterate yo was chosen to be equal to the 
zeroth WKB approximation ik, one can state that the p- 
th QLM iterate contains 2^ exact terms. In addition the 
expansion of each QLM iterate has the proper structure 
whose terms are identical to those of the WKB series but 
have only approximately correct coefficients. 

The 2P law is, of course, not accidental. The QLM 
iterates are quadratically convergent [l|, Q, that is the 
norm of the difference of the exact solution and the p-th 
QLM iterate Ijj/ — 2/p|| is proportional to the square of 
the norm of the differences of the exact solution and the 
{p - l)-th QLM iterate: 



\\y-yp\\ - \\y-yp-i\y 



(6) 



Here the norm \\g\\ of the function g{x) is the maximum 
of the function g{x) on the whole interval of values of x. 
Since j/o contains just one correct WKB term of power g^ 
and thus Hy — 2/o|| is proportional to g, one has to expect 
that 1 1 y— 2/1 1 1 ^ g'^ and thus j/i contains two correct WKB 
terms of powers g'^ and g^. The difference \\y — y2\\ ^ 
\\y ~ yi\\ ^ 5^ so that y2 contains four correct WKB 
terms of powers g", g^, g^ and g^ . Finally, the difference 
\\y " J/3 1 1 should be proportional to g^, and therefore 2/3 
has to contain eight correct terms with powers between 
g° and g"^ inclusively. This explains the 2^ law. 

The exact quantum mechanical quantization condition 
for the energy |3,ll3 '^^s the form: 



J 



y{z)dz — i 27rn. 



(7) 



Here y{z) is the logarithmic derivative of the wave func- 
tion, given by Eq. Q, z = gr, n — 0,1, 2, ... is the bound 
state number which counts the number of poles of y{z); 
the integration is along a path C in the complex plane 
encircling the segment of the 5Rz axis between the turning 
points. 

The p-th QLM iterate yp{z), as we have seen, con- 
tains, in addition to 2^ exact WKB terms of powers 17°, 
g^ , . . ., g^ ~^, also an infinite number of structurally cor- 
rect WKB terms of higher powers of g with approximate 
coefficients. One can expect therefore that the quantiza- 
tion condition with y{z) approximated by any QLM 
iterate yp{z), p = 1,2, . . ., including the first, gives more 
accurate energy than the usual WKB quantization con- 
dition which is obtained by substituting into exact quan- 
tization condition |(7} the WKB expansion up to the first 
power ol g r^ h, that is y{z) = ik{z) — {dk{z)/dz)/{2k{z)), 
and neglecting all higher powers of g in the expansion. 
Indeed, we will prove now that Eq. |(7} with y{z) ap- 
proximated by any QLM iterate leads to exact energies 



not only for the Coulomb and harmonic oscillator poten- 
tials as it was shown earlier in Ref. [6j|, but for many 
other well known physical potentials used in molecular 
and nuclear physics such as the Poschl- Teller, Hulthcn, 
Hylleraas, Morse, Eckart etc. The WKB quantization 
condition yields the exact energy only for the first two 
potentials, but not for the rest of them. 

Let us prove it on the example of the Hulthen poten- 
tial V{r) = -Aexp{-r/a)/{l - exp {-r/a)),A > 0, < 
r < oo, which plays an important role in molecular and 
nuclear physics. 

To compute energy levels in the quasilinear approxi- 
mation we have to use the QLM equation (0) which after 
switching to the variable t — exp(— r/a),0 < t < 1 has 
the form 



t dyp 
a dt 



— -;^ = yl-'^yp vp-i + f- A 



t 



1 -t 



(8) 



Here e — ~E and E is the energy. For convenience of 
further computations we set here 2mjf? equal to unity. 
The quantization condition Q in variable t is given by 
Jp — a§ fj {yjy{t)/t) dt — i2Trn, p ~ 1,2, . . . 

At the singular point t ~ of the integrand, Eq. 
reduces to —{t/a){dyp/dt) — y^ — 2yp yp-i + e whose so- 
lution is j/p — Cp, where Cp is a constant satisfying the 
algebraic equation Cp — 2cp Cp_i -I- e = 0. Since at large p 
we expect j/p-i —> yp ^ y, where y is an exact solution 
at i = 0, it means, in view of Cp being a constant, that 
we should have Cp = Cp_i = c, that is, we are looking for 
a fixed point solution of this algebraic equation, which is 
Cp = y/e. The positive sign before the root is chosen since 
the first term in expansion of yp{t) in the WKB terms is 
ik{t). Thus yp{t) 
2/p(0) ~ +,re. 



ik{t) = i^/-e+ At/{1 -t), that is 

In the same way one finds that the residues at other 
singular points t — 1 and i = oo equal 1/a and Ve + A, 
respectively. After taking into account that poles t ^ 
and t = 1 are encircled in the opposite direction com- 
pared with the pole at infinity (which is encircled coun- 
terclockwise) and the reinstatement of the factor 21X1 jf? 
the exact quantization condition gives yj2/ma?/W{'-- %/e+ 
i/e -t- A) = n + 1, 71 = 0, 1, 2, . . .. This expression coin- 
cides with the exact one for the Hulthen potential calcu- 
lated in [ill n^ and is different from the WKB quanti- 
zation condition yj 2ma? / tt^ {— ^ + \/e + A) = n + ^. 

Similar computations show that the quantization con- 
dition lO in any QLM approximation including the first 
leads to exact energies for all the potentials mentioned 
above and for other potentials with a simple singular 
structure. 

For more complicated potentials numerical calculation 
is necessary. However, as we will see now, already the 
first QLM iterate, given by the closed analytic expression, 
is extremely accurate. For the zeroth iterate yo{z) one 
can choose the usual WKB approximation. However, this 
choice has unphysical turning point singularities. Conse- 
quently, if yp{z) in Eq. |(2Jl is a discontinuous function of 
z in a certain interval, then [2| yp+i{z) or its derivatives 



could also be discontinuous functions in this interval, so 
the turning point singularities oiyo{z) will unfortunately 
propagate to the next iterate. To avoid this we choose 
the Langer WKB wave function 13] as the zeroth iter- 
ation. This function near the turning points a and b is 
given by the simple analytic expression 






jmM 




dSfir) 



\k{s)\ds 



(9) 



Here Ai denotes the Airy function, i = a,b, k{r) = 
{2m/h'^)[{E-V{r))-{l + l/2)^/{2mr^)], d is -1 for 
a < r < b and 1 for r < a,r > b, and Ca — I, Cb — (—1)", 
where n — 0, 1, 2, . . . is the number of the bound state. 
It is easy to check that Xa{.r) and Xb{f) coincide at some 
point at the interval (a, b) between the turning points, 
are continuous across them and coincide with the usual 
WKB solutions far from them. 

Let us consider a couple of examples. The exact en- 
ergy of the ground state of the anharmonic oscillator 
V{r) = r'^ is 2.393 644 016 482 303 115 6 in atomic units 
with mass set to unity, m = 1. This result is obtained by 
us using the Runge-Kutta method and quadruple preci- 
sion arithmetic. The WKB energy is different by 2.8% 
and equals 2.32662, while the first-iteration QLM energy 
equals 2.39475 and differs from the exact energy only by 
0.046%. The QLM energy coincides with the exact en- 
ergy in all twenty digits after the sixth iteration. 

The graph in Fig.^displays the logarithm of the differ- 
ence between the exact and WKB solutions and between 
the exact solution and the first QLM iteration. One can 
see that the difference between the exact solution and the 
first QLM iteration is two orders of magnitude smaller 
than the difference between the exact and the WKB so- 
lutions, that is, just one QLM iteration increases the ac- 
curacy of the wave function by two orders of magnitude 

The second example is the modified Coulomb potential 



V{r) = -- 
2p 



l{l 



1) - ia2 



4" 



p2(p- 



aEr 



obtained when the equal-mass two-body Dirac equation 
with the static Coulomb interaction is reduced to the 
Schrodinger equation [lj|. The exact energy of the ^Sq 
ground state is 0.999 993 340 148 538 880 12 in atomic 
units with double mass set to unity, 2M — 1. This re- 
sult was obtained in the work [l5j by an elaborate com- 
putation using the finite element method and verified 
by ourselves using the Runge-Kutta method in quadru- 
ple precision. WKB in this case predicts the energy 
very accurately since the potential is very close to the 
Coulomb potential, for which WKB energy is precise. 
It equals 0.999 986 680 and differs from the exact one 
only by 6.66 x 10~^. The first-iteration QLM energy 



log;olZof-»:J-XexactWi 




FIG. 1: Logarithm of the difference between the exact so- 
lution Xcxact and the WKB solution xo (dashed curve) and 
between Xoxact and the first QLM iterate Xm„ (solid curve) 
for the ground state of the anharmonic oscillator [lq|. 



equals 0.999 993 335 4 and differs from the exact energy 
by 5 X 10~^, that is, it is more accurate than the WKB 
energy by three orders of magnitude. The QLM energy 
coincides with the exact one up to the twentieth digit 
after the sixth iteration. 

The calculation shows that the difference between the 
exact wave function and the first QLM iteration is, as 
in previous example, by two orders of magnitude smaller 
than the difference between the exact and the WKB solu- 
tions. Thus also in this case one QLM iteration increases 
the accuracy of the wave function by a remarkable two 
orders of magnitude. 

In conclusion, we have shown that the quasilineariza- 
tion method (QLM) which approaches solution of the 
Riccati-Schrodinger equation by approximating the non- 
linear terms by a sequence of the linear ones, and is not 



based on the existence of a small parameter, sums the 
WKB series. The expansion of the p-th QLM iteration 
in powers of h reproduces the structure of the WKB series 
generating an infinite number of the WKB term with 2^ 
terms of the expansion reproduced exactly and a similar 
number approximately. As a result one expects that the 
exact quantization condition with integrand replaced by 
any QLM iterate including the first gives more accurate 
energy than the WKB quantization condition which is 
obtained by substituting into Eq. Q of the WKB expan- 
sion up to the first power of h and neglecting all higher 
powers of H. We show on examples of the Hulthen po- 
tential that it in fact, the QLM energy is exact already 
in the first iteration. Similarly, one can show that the 
approximation by the first QLM iterate in Eq. (71) leads 
to exact energies for many well known physical potentials 
such as the Coulomb, harmonic oscillator, Poschl- Teller, 
Hulthen, Hylleraas, Morse, Eckart etc. For other poten- 
tials which have more complicated analytical structure 
we show on examples of the anharmonic oscillator and 
modified Coulomb potentials that the use of the Langer 
WKB solution as an initial guess already in the first QLM 
approximation gives energies and wave functions at least 
two orders of magnitude more accurate than the WKB 
results. Such a QLM solution, unlike the usual WKB so- 
lution, displays no unphysical turning point singularities. 
Since the first QLM iterate is given by a close analytic 
expression it allows one to estimate analytically the role 
of different parameters and their infiuence on properties 
of a quantum system with much higher precision than 
provided by the WKB approximation. In addition, it 
was shown that six QLM iterations are usually enough 
to obtain both the wave function and the energy with 
extreme accuracy of twenty significant figures. 
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The dips are artifacts of the logarithmic scale, as the log- 
arithm of the absolute value of the difference of two solu- 
tions goes to —CO at points where the difference changes 
sign. The overall accuracy of the solution can therefore 
be inferred only at x values not too close to the dips. 



